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Abstract. Globular clusters offer ideal laboratories to test the predictions of stellar evolution. 
When doing so with spectroscopic analyses during the 1990s, however, the parameters we de- 
rived for hot horizontal branch stars deviated systematically from theoretical predictions. The 
parameters of cooler, A-type horizontal branch stars, on the other hand, were consistent with 
\l ' evolutionary theories. In 1999, two groups independently suggested that diffusion effects might 
, cause these deviations, which we verified subsequently. I will discuss these observations and 
analyses and their consequences for interpreting observations of hot horizontal branch stars. 
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1. Historical Background 

Globular clusters are densely packed, gravitationally bound systems of several thou- 
sands to about one million stars. The dimensions of the globular clusters are small com- 
, pared to their distance from us: half of the light is generally emitted within a radius of 
' less than 10 pc, whereas the closest globular cluster has a distance of 2 kpc and 90% 
lie more than 5 kpc away. We can thus safely assume that all stars within a globular 
cluster lie at the same distance from us. With ages in the order of 10^° years globular 
, clusters are among the oldest objects in our Galaxy. As they formed stars only once in 
the beginning and the duration of that star formation episode is short compared to the 
current age of the globular clusters the stars within one globular cluster are essentially 
coeval. In addition all stars within one globular cluster (with few exceptions) show the 
same initial abundance pattern (which may differ from one cluster to another). Globular 
clusters are thus the closest approximation to a physicist's laboratory in astronomy. 



, The horizontal branch, which is the topic of this article, was discovered by ten Bruggencate (1927)[ 



when he used Shapley's (1915) data on M 3 and other clusters to plot magnitude 
versus colour (replacing luminosity and spectral type in the Hertzsprung-Russell dia- 
gram) and thus produced the first colour-magnitude diagrams ( "Farbenhelligkeits- 
diagramme" ) . In these colour-magnitude diagrams (CMD's) ten Bruggencate noted the 
presence of a red giant branch that became bluer towards fainter magnitudes, in agree- 



ment with Shapley (1915) In addition, however, he saw a horizontal branch ("hori- 
ZONTALER Ast" ) that parted from the red giant branch and extended far to the blue 
at constant brightness. As more GMD's of globular clusters were obtained it became 
obvious that the relative numbers of red and blue horizontal branch stars (i.e. the hor- 
izontal branch morphology) varied quite considerably between individual clusters, with 
some clusters showing extensions of the blue horizontal branch (so-called "blue tails") 
towards bluer colours and fainter visual magnitudes, i.e. towards hotter temperaturesf. 

t The change in slope of the horizontal branch towards higher temperatures is caused by the 
decreasing sensitivity of _B — V to temperature on one hand and by the increasing bolometric 
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Sandage & Wallerstein (1960) noted a correlation between the metal abundance and the 



horizontal branch morphology seen in globular cluster CMD's: the horizontal branch 
(HB) became bluer with decreasing metallicity. 



About 25 years after the discovery of the horizontal branch Hoyle & Schwarzschild (1955) 



were the first to identify horizontal branch stars with post-red giant branch stars that 
burn helium in the central regions of their cores. fFaulkner (1966) managed for the first 
time to compute zero age horizontal branch (ZAHB) models that qualitatively re- 
produced the observed trend of HB morphology with metallicity without taking into 
account any mass loss but assuming a rather high helium abundance of Y = 0.35. 



Iben & Rood (1970) however, found that they could ". . . account for the observed spread 
in colour along the horizontal branch by accepting that there is also a spread in stellar 
mass along this branch, bluer stars being less massive ( on the average ) and less lumi- 
nous than redder stars." Comparing HB models to observed globular cluster CMD's 



Rood (1973) found that an HB that . . is made up of stars with the same core mass 



and slightly varying total mass, produces theoretical c-m diagrams very similar to those 
observed. ...A mass loss of perhaps 0.2 Mq with a random dispersion of several hun- 
dredths of a solar mass is required somewhere along the giant branch. " The assumption of 
mass loss on the red giant branch diminished the need for very high helium abundances. 

Thus our current understanding sees horizontal branch stars as stars that burn helium 
in a core of about 0.5 and hydrogen in a shell and evolve to the asymptotic giant 
branch, when the helium in the core is exhausted (for a review on HB evolution see 
Sweigart 1994). The more massive the hydrogen envelope is the cooler is the resulting 
star at a given metallicityf. The masses of the hydrogen envelopes vary from ^0.02 
at the hot end of the horizontal branch (about 30,000 K) to 0.3-0.4 M© for the cool HB 
stars at about 4000-5000 K (depending on metallicity, e.g. Dorman et al. 1993). The 
stable red and blue HB stars are separated by the variable RR Lyrae range at about 
6500-7500 K. This article deals with blue HB stars, which at effective temperatures of 
about 8,000 K to 20,000 K show spectra rather similar (at moderate resolution) to main 
sequence stars of spectral types A and B and are therefore called HBA and HBB stars. 
In the field of the Milky Way such stars are often denominated by FHB (field HB star) 
and used as tracers for halo structure. 



2. Atmospheric parameters 

Already early studies of HBA and HBB stars in globular clusters showed discrepancies 



between observational results and theoretical expectations: Graham & Doremus (1966) 
mentioned that the comparison of (ci)o vs. (6 — y)o for 50 blue HB stars in NGC 6397 to 
models from Mihalas (1966) indicated low surface gravities and a mean mass of O.3M0 



(O.4M0 ) for solar (neghgible) hehum abundance, assuming (to — M)o = 12.0 and Eb-v 
= 0.16. Later spectroscopic analyses of HB stars (see cited papers for details) in globular 
clusters reproduced this effect (cf. Fig.P): Crocker et al. (1988, M 3, M 5, M 15, M 92, 
NGC 288), de Boer et al. (1995, NGC 6397), Moehler et al. (1995, 1997a, M 15), Moehler 
et al. (1997b, NGC 6752). 

The zero age HB (ZAHB) in Fig.^marks the position where the HB stars have settled 
down and started to quietly burn helium in their cores. The terminal age HB (TAHB) is 
defined by helium exhaustion in the core of the HB star {Yc < 0.0001). For temperatures 

correction for hotter stars (i.e. the maximum of stellar flux is radiated at ever shorter wavelengths 
for increasing temperatures, making stars fainter at V) on the other hand. 

f Due to the higher opacities in their envelopes metal-rich HB stars are cooler than metal-poor 
ones with the same envelope mass. 
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Figure 1. The results of Crocker et al. (1988, M 3, M 5, M 92, NGC 288), de Boer et al. 
(1995, NGC 6397), Moehler et al. (19 95, 1997a, M 15), and Moehler et al. (1997b, NGC 6752) 
compared to evolutionary tracks from [Dorman et al. (1993)] ZAHB and TAHB stand for zero 
age and terminal age HB (see text for details). I'he short-dashed lines mark the regions of low 
logg (see text for details). 



between 12,000 K and 20,000 K the observed positions in the (logg. Toff )-diagram fall 
mostly above the ZAHB and in some cases even above the TAHB. This agrees with the 
finding of Saffer et al. (1997) that field HBB stars show a larger scatter away from the 
ZAHB in Teff, logg than the hotter subdwarf B stars with Toff > 20,000 K. Knowing the 
atmospheric parameters of the stars and the distances to the globular clusters we can 
determine masses for the stars (cf. Moehler et al. 1995, 1997b, de Boer et al. 1995). While 
the stars in M 3, M 5, and NGC 6752 have mean masses consistent with the canonical 
values, the hot HB stars in all other clusters show masses that are significantly lower 
than predicted by canonical HB evolution ~ even for temperatures cooler than 12,000 K. 

Also some UV observations suggest discrepancies between theoretical expectations and 
observational results: The lUE (International Ultraviolet Explorer) and HUT (Hopkins 
Uhraviolet Telescope) spectra of M 79 (Ahner & Matilsky 1993, Dixon et al. 1996) 
suggest lower than expected gravities and higher than expected metallicities for hot HB 
stars (but see Vink et al. 1999, who do not need low surface gravities to fit the HUT 
data). I Hill et al. ( 19967 find fr om Ultraviolet Imaging Telescope (UIT) photometry of 
M 79 that stars bluer than mi52 — 771249 = —0.2 lie above the ZAHB, whereas cooler 
stars scatter around the ZAHB. UIT data of M 13 (Parise et al. 1998) find a lack of stars 
close to the ZAHB at a colour (temperature) range similar to the low logg range shown 
in Fig. n] Landsman et al. (1996) on the other hand find good agreement between UIT 
photometry of blue stars in NGC 6752 and a standard ZAHB. 
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3. Atmospheric abundances 

It has been realized early on that the blue HB and blue tail stars in globular clus- 
ters show weaker helium lines than field main sequence B stars of similar temperatures: 
Searle & Rodgers (1966, NGC 6397); Greenstein fc Miinch (1966, M 5, M 13, M 92); 



Sargent (1967, M 13, M 15, M 92). Greenstein et al. (1967) aheady suggested diffu 



sion to explain this He deficiency. Michaud et al. (1983) performed the first theoretical 



study of diffusion effects in hot horizontal branch stars. Using the evolutionary tracks of 



Sweigart & Gross (1976) they found for the metal-poor models that "in most of each en- 
velope, the radiative acceleration on all elements (i.e. C, N, O, Ca, Fe) is much larger than 
gravity which is not the case in main-sequence stars. " The elements are thus pushed to- 
wards the surface of the star. Turbulence affects the different elements to varying extent, 
but generally reduces the overabundancesf. Models without turbulence and/or mass loss 
(which may reduce the effects of diffusion) predict stronger He depletions than observed. 
A weak stellar wind could alleviate this discrepancy (Heber 1986, Michaud et al. 1989, 
Fontaine & Chayer 1997, and Unglaub & Bues 1998 discuss this effect, albeit for hotter 
stars). The extent of the predicted abundance variations varies with effective tempera- 
ture, from none for HB stars cooler than about 5800±500K (due to the very long diffusion 
timescales) to 2 - 4 dex in the hotter stars (the hottest model has Tcff — 20,700 K) and 
also depends on the element considered. 

Observations of blue HB and blue tail stars in globular clusters support the idea of 
diffusion being active above a certain temperature: Abundance analyses of blue HB stars 
cooler than 11,000 K to 12,000 K in general show no deviations from the globular cluster 
abundances derived from red giants, while for hotter stars departures from the general 
globular cluster abundances are found, e.g. iron enrichment to solar or even super-solar 
values and strong helium depletion (see Moehler 2001 for references and more details). 



This agrees with the finding of .Altner fc Matilsky (1993) and Vink et al. (1999) that 



solar metallicity model atmospheres are required to fit the UV spectra of M 79. 

All this evidence supports the suggestion of [Grundahl et al. (1999) that the onset of 



diffusion in stellar atmospheres may play a role in explaining the jump along the HB 
towards brighter u magnitudes at effective temperatures of about 11,500 K. This jump 
in u,w — y is seen in all CMD's of globular clusters that have Stromgren photometry 
of sufficient quality. The effective temperature of the jump is roughly the same for all 
clusters, irrespective of metallicity, central density, concentration or mixing evidence, 
and coincides roughly with the onset of the "low logg problem" seen in Fig. ^ at ToS 
«11,000 K to 12,000 K. This in turn coincides with the region where surface convection 
zones due to hydrogen and Hei ionization disappear in HB stars (Sweigart 2002). 

Radiative levitation of heavy elements decreases the far-UV ffux and by backwarm- 
ing increases the ffux in u. Grundahl et al. (1999) show that the use of metal-rich at- 
mospheres ([Fe/H] = +0.5 for scaled-solar ATLAS9 Kurucz model atmospheres with 
logEFe,© = 7.60) improves the agreement between observed data and theoretical ZAHB 
in the u,u — y-CMD at effective temperatures between 11,500 K and 20,000 K, but it 
worsens the agreement between theory and observation for hotter stars in the Stromgren 
CMD of NGC 6752 (see their Fig. 8). Thus diffusion may either not be as important 
in the hotter stars or the effects may be diminished by a weak stellar wind. The gap 



t [Micha ud (1982)1 and [Charb onncau & Michaud (1988) showed that meridional circulation 
can prevent gravitatijmal settling and that the limiting rotational velocity decreases with de- 
creasing logg. Bchr et al. (2000b) note that two of the HB stars hotter than 10,000 K show 
higher rotational velocities and much smaller abundance deviations than other stars of similar 
temperature. 
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at {B — V)o ~ discussed by Caloi (1999) is not directly related to the u-jump as it 
corresponds to an effective temperature of about 9000 K and is also not seen in every 
cluster (which would be expected if it were due to an atmospheric phenomenon) . The gap 
at Teff « 13,000 K seen in the ci,b — y diagram of field horizontal branch stars (Newell 
1973, Newell & Graham 1976) may be related to the u-jump as the Ci index contains u. 

The abundance distribution within a stellar atmosphere influences the temperature 
stratification and thereby the line profiles and the flux distribution of the emergent spec- 
trum. A deviation in atmospheric abundances of HB stars from the cluster metallicity 
due to diffusion would thus affect their line profiles and fiux distribution. Model atmo- 
spheres calculated for the cluster metallicity may then yield wrong results for effective 
temperatures and surface gravities when compared to observed spectra of HB stars. This 
effect could explain at least part of the observed discrepancies (see Sect. for more de- 
tails). Self-consistent model atmospheres taking into account the effects of gravitational 
settling and radiative levitation are, however, quite costly in CPU time and have started 
to appear only quite recently for hot stars (Dreizler & Wolff 1999, Hui-Bon-Hoa et al. 
2000). 



4. Rotational velocities 

Peterson (1983, 1985a, 1985b) found from high-resolution spectroscopic studies that 
clusters with bluer HB morphologies show higher rotation velocities among their HB 
stars. However, the analysis of |Peterson et al. (1995) shows that while the stars in M 13 
(which has a long blue tail) rotate on average faster than those in M 3 (which has only a 
short blue HB), the stars in NGC 288 and M 13 show slower rotation velocities at higher 
temperatures. These results are consistent with those reported by Behr et al. (2000ajl 
who determined rotational velocities for stars as hot as 19,000 K in M 13. They found 
that stars hotter than about 11,000 K have significantly lower rotational velocities than 
cooler stars and that the change in mean rotational velocity may coincide with the gap 
seen along the blue HB of M 13 (see Ferraro et al. 1998 for an extensive discussion of 
gaps). Also the resuffs of Cohen & McCarthy (1997, M 92) and Behr et al. (2000b, M15) 
show that HB stars cooler than ^11,000 K to 12,000 K in general rotate faster than 
hotter stars. 

The studies by Behr (2003a) and Recio et al. (2002, 2004), which both consider several 
clusters, confirm the trends mentioned above. Below the diffusion threshold about 20% 
to 30% of the blue HB stars show rotation velocities of vs'mi 20 ... 30 km/s, whereas 
there are no fast rotators among the hotter stars. From observations of field HB stars 



Behr (2003b) and Carney et al. (2003) note that red and cool blue HB stars show similar 



distributions of rotational velocities (after accounting for the larger radii of the red HB 
stars), whereas field RR Lyrae stars show no evidence for rotation. 



Sills & Pinsonneault (2000) study theoretical models for the rotation of HB stars and 



find that the observed rotation of cool blue HB stars in M 13 can be explained if their red 
giant precursors have rapidly rotating cores and differential rotation in their convective 
envelopes and if angular momentum is redistributed from the rapidly rotating core to 
the envelope on the horizontal branch. If, however, turn-off stars rotate with less than 
4 km/s, a rapidly rotating core in the main-sequence stars (violating helioseismological 
results for the Sun) or an additional source of angular momentum on the red giant branch 
(e.g. mass transfer in close binaries or capture of a planet as described by Soker & Harpaz 



2000) are required to explain the rotation of blue HB stars. Talon & Charbonnel (2004) 



speculate that such differential rotation could be understood within the framework of 
internal gravitational waves. The change in rotation rates towards higher temperatures is 
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not predicted by the models of Sills & Pinsonneault (2000) but could be understood as a 
result of gravitational settling of helium, which creates a mean molecular weight gradient, 
that then inhibits angular momentum transport to the surface of the star. Sweigart (2002) 



suggests that the weak stellar wind invoked to reconcile observed abundances in hot HB 
stars with diffusion calculations (cf. Sect. ^ could also carry away angular momentum 
from the surface layers and thus reduce the rotational velocities of these stars. 



Soker & Harpaz (2000) argue that the distribution of rotational velocities along the 



HB can be explained by spin-up of the progenitors due to interaction with low-mass 
companions, predominantly gas-giant planets, in some cases also brown dwarfs or low- 
mass main-sequence stars (esp. for the very hot HB stars). The slower rotation of the 
hotter stars in their scenario is explained by mass loss on the HB, which is accompanied 
by efficient angular momentum loss. This scenario, however, does not explain the sud- 
den change in rotational velocities and the coincidence of this change with the onset of 
radiative levitation. 



5. Where do we stand? 

Analysis of a larger sample of hot HB stars in NGC 6752 (Moehler et al. 2000) showed 
that the use of model atmospheres with solar or super-solar abundances removes much 
of the deviation from canonical tracks in effective temperature, surface gravity and mass 
for hot HB stars discussed in Sect. |21 However, some discrepancies remain, indicating 
that the low log g, low mass problem cannot be completely solved by scaled-solar metal- 
rich atmospheres (which do reproduce the u-jump reported by Grundahl et al. 1999). As 



Michaud et al. (1983) noted diffusion will not necessarily enhance all heavy elements by 



the same amount and the effects of diffusion vary with effective temperature. Elements 
that were originally very rare may be enhanced even stronger than iron. The question of 
whether diffusion can fully explain the "low gravity" problem cannot be answered without 
detailed abundance analyses to determine the actual abundances and subsequent analy- 
ses using model atmospheres with non-scaled solar abundances (like ATLAS12, Kurucz 



1992). Additional caution is also recommended by the results of Moehler et al. (2003) 



on M 13, where even the use of the metal-rich model atmospheres does not eliminate 
the problem of the low masses or the low gravities for effective temperatures between 
12,000 K and 16,000 K. In that cluster we also found significant disagreement between 
atmospheric parameters derived from Stromgren photometry or from Balmer line fitting 
for stars cooler than about 9,000 K, for which we found no explanation. 

Still unexplained are also the low masses found for cool blue HB stars (which are not 
affected by diffusion) in, e.g., NGC 6397 and M 92. For those stars a longer distance 
scale to globular clusters would reduce the discrepancies (Moehler 1999). Such a longer 
distance scale has been suggested by several authors using Hipparcos results for metal- 
poor field subdwarfs to determine the distances to globular clusters by fitting their main 
sequence with the local subdwarfs. One should note, however, that de Boer et al. (1997)[ 
report that Hipparcos parallaxes for field HBA stars still yield masses significantly 
below the canonical mass expected for these objects. Carretta et al. (2000) present an 
extensive and excellent discussion of various globular cluster distance determinations and 
the zoo of biases that affect them, while Gratton et al. (2003) concentrate on the error 
budget of distances from main sequence fitting in their paper. 

The problem of the different rotation velocities for cool and hot HB stars, however, 
remains rather stubbornly opposed to any attempted solution. 
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Discussion 

Freytag: What are "real" abundances in contrast to "diffusion" abundances? Can they 
still be reliably defined? 

MoEHLER: By real abundances I mean the abundances, which the star would show if its 

atmosphere were not affected by diffusion (like the abundances derived from red giants 
in globular clusters). And those cannot be derived if the spectrum is affected by diffusion 

Preston: Have you compared the abundances of hot HB stars with those derived from 
red giants, which are affected in a different way by diffusion processes? 

MOEHLER: Red giants should show no diffusion effects in their atmospheres due to their 
deep convection zones. HB stars with diffusion usually show strong enhancements of 
heavy elements compared to red giants in the same cluster. 

CORBALLY: Some dozen years ago I was working with Richard Gray on field horizontal 
branch stars identified by Dave Philip. These were all slightly cooler than those you have 
shown, i.e. truly A-type horizontal branch stars. We found that they lay near or slightly 
above the main sequence, indeed they were metal-poor as you said, and that two, from 
their Balmer line profiles at medium resolution, were probably helium-rich. Could you 
comment on the last possibility? 

MOEHLER: Not really, as I find this a very puzzling result, which cannot be easily un- 
derstood in the terms of diffusion. 



